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Both the gauge- invariant fermion Green function and gauge-dependent conventional Green func- 
tion in 2 + 1 dimensional QED are studied in the large TV limit. In temporal gauge, the infra-red 
divergence of gauge-dependent Green function is found to be regulariable, the anomalous dimen- 
sion is found to be ri = . ®f . This anomalous dimension was argued to be the same as that of 
gauge-invariant Green function. However, in Coulomb gauge, the infra-red divergence of the gauge- 
dependent Green function is found to be un-regulariable, anomalous dimension is even not defined, 
but the infra-red divergence is shown to be cancelled in any gauge-invariant physical quantities. The 
gauge-invariant Green function is also studied directly in Lorentz covariant gauge and the anomalous 
dimension is found to be the same as that calculated in temporal gauge. 



I. INTRODUCTION 

In high energy physics, 2 + 1 dimensional massless 
Quantum Electro-Dynamics ( QED3 ) was studied as an 
interesting theoretical toy model to gain insight into dy- 
namical symmetry breaking of 3 + 1 dimension QED [1] . 
In condensed matter systems, QED3 may appear as an 
effective low energy theory describing doped Mott insula- 
tors [2-4] . The connection of dynamical mass generation 
in QED3 with Anti-ferromagnets was also discussed in 
[3] . QED3 with possible Chern-Simon term may also de- 
scribe Fractional Quantum Hall (FQH) transitions [5-9] . 

In the two component notations suitable for describ- 
ing FQH transitions [6-9], the standard 2 + 1 dimen- 
sional massless Quantum Electro-Dynamics ( QED3 ) 
Lagrangian in Euclidean space is: 

£ = ^o7m(^ _ iea^ipa + ^(/^) 2 (!) 

Where the Dirac matrices in three dimensions are 7 M = 
ia^ with <jy.{\i — 1,2,3) being three Pauli matrices [7] 
and a — 1, • • • , N are N species of Dirac fermion. 

The conventional Green function in QED3 is defined 
as: 

G(x 1 ,x 2 )=<^(x 1 )tp(x 2 )> (2) 

In momentum and real space, the fully-interacting 
Green function in Eqn.2 takes the form: 

G(*) = ^gf, G(*) = 3gf (3) 

where the convention we chose for the sign of the anoma- 
lous dimension 77 is that a positive rj slows down the al- 
gebraic decay of the Green function in real space, 77 can 
be calculated by standard Renormalization Group (RG) 
by extracting UV divergences. Unfortunately, this Green 
function is not gauge invariant, 77 depends on the fixed 
gauge in which the calculation is done. 

The gauge invariant Schwinger Green function which 
is the focus of this paper is 



G™( Xl ,x 2 ) =< Tp( Xl )e ie ^ aAm ^(x 2 ) > (4) 

where the inserted Dirac string makes the Schwinger 
Green function gauge invariant. 

Qinv depends on the the integral path C from x\ to x 2 . 
For simplicity reason, we take C to be a straight line. 

As far as I know, the importance of gauge-invariant 
Green function in condensed matter system was first real- 
ized in FQH system [11]. By a singular gauge transforma- 
tion which attaches two flux quantums to each electron, 
an electron in an external magnetic field was mapped 
to a composite fermion in a reduced magnetic field. Al- 
though transport properties which are directly related to 
two particle Green functions can be directly studied in 
the composite fermion language, the tunneling density of 
states which is directly related to the single particle elec- 
tron Green function is much more difficult to study. In 
fact, the single particle electron Green function is equal 
to the gauge-invariant Green function of the composite 
fermion which was evaluated for non-relativistic fermions 
by phenomelogical arguments in Ref. [11]. 

Most recently, the importance of gauge-invariant 
Green function of fermion to Angle Resolved Photo- 
Emission (ARPES) data in high temperature supercon- 
ductors was discovered independently by Rantner and 
Wen (RW) [12] and the author [14] in different contexts. 
Starting from SU(2) gauge theory of doped Mott in- 
sulators [4], RW discussed the relevance of this gauge- 
invariant Green function to ARPES data. They also 
pointed out that in temporal gauge, the equal-space 
gauge-invariant Green function in Eqn.4 is equal to that 
of the conventional gauge dependent one in Eqn.2. Fur- 
thermore, they calculated the anomalous dimension of 
the conventional Green function in temporal gauge by 
large N expansion. However, both the sign and mag- 
nitude of their result differs from that achieved in this 
paper. Starting from a complementary (or dual ) ap- 
proach pioneered by Balents et ed [13], the author studied 
how quantum [15] or thermal [14] fluctuations generated 
hc/2e vortices can destroy d-wave superconductivity and 
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evolve the system into underdoped regime at T = or 
pseudo-gap regime at finite T. The author found that 
only in weakly type II cZ-wave superconductors where the 
interaction between vortices is short-ranged, the effective 
low energy theory is described by QED3. For underdoped 
cuprates which are between quantum anti-ferromagnet 
and strongly type II d wave superconductors, the correct 
T = theory is still unknown. However, in the vortex 
plasma regime around the finite temperature Kosterlizc- 
Thouless transition [26], the vortices can be treated by 
classical hydrodynamics. By Anderson singular gauge 
transformation which attaches the flux from the classi- 
cal vortex to the quasi-particles of (Z-wave superconduc- 
tors [16-19], the quasi-particles (spinons) are found to 
move in a random magnetic field generated by the clas- 
sical vortex plasma [14]. The electron spectral function 
G(x, t) —< C Q (0, 0)C^(x, t) > is the product of the clas- 
sical vortex correlation function and the gauge invariant 
Green function of the spinon in the random magnetic 
field. Obviously, this static gauge invariant Green func- 
tion is different from that calculated in this paper and 
will be investigated by completely different methods in 
separate publication [20]. 

In this paper, for pure theoretical interests, but with 
no application to high temperature superconductors in 
mind, we study the gauge-invariant fermion Green func- 
tions Eqn.4 in large N limit. By applying the methods 
developed to study clean [8] and disordered [9] FQH tran- 
sitions and superconductors to insulator transitions [10], 
we study the gauge-dependent Green function both in 
temporal gauge and in Coulomb gauge. In both gauges, 
we regularize the UV divergences in two different cut-offs. 
In both gauges, we also run into expected IR singular- 
ities. We also make detailed comparisons between the 
natures of the IR singularities in these two gauges. In 
temporal gauge, the infrared divergence is in the mid- 
dle of the contour integral along the real axis. From 
physical prescription, it can be regularized by deforming 
the contour into complex plane, the anomalous dimen- 
sion is found to be n = ■ From RW's observation 
that equal-space gauge invariant Green function is the 
same as the equal space gauge dependent one in tempo- 
ral gauge, we expect this anomalous dimension to be the 
same as that of gauge-invariant Green function. How- 
ever, in Coulomb gauge, the infra-red divergence is at 
the two ends of the contour integral along the real axis, 
therefore un-regulariable, anomalous dimension is even 
not defined. This should not cause too much concern, be- 
cause the anomalous dimension in Coulomb gauge does 
not corresponds to any physical gauge-invariant quan- 
tity. This infra-red divergence can be shown to be can- 
celed in any physical gauge invariant quantities such as (3 
function and correlation length exponent v as observed 
in Refs. [8]. We also study the gauge- invariant Green 
function directly by Lorentz covariant calculation with 
different gauge-fixing parameters and find that the expo- 



nent is independent of the gauge fixing parameters and 
is exactly the same as that found in temporal gauge. 

The paper is organized as follows. In the next section, 
we briefly review the previous work on gauge dependent 
Green function in Lorentz covariant Landau gauge [1]. 
In section III, we calculate the gauge dependent Green 
function in temporal gauge with two different cut-offs. In 
section IV, we discuss the gauge dependent Green func- 
tion in Coulomb gauge also with two different cut-offs. 
We make detailed comparisons with the calculations done 
in temporal gauge. In section V, we calculate directly 
the gauge invariant Green function in Lorentz covariant 
gauge with different gauge fixing parameters. Finally, we 
reach conclusions. 



II. THE CALCULATION IN LANDAU GAUGE 

In this section, we will briefly review the previous work 
on the gauge-dependent Green Eqn.2 in QED3. There is 
nothing new in this section. But it is constructive to com- 
pare the calculations done in Lorentz covariant gauge in 
this section with those done in temporal gauge in section 
III and in Coulomb gauge in section IV. It is also a useful 
first step for the Lorentz covariant calculation of gauge- 
invariant Green function Eqn.4 to be presented in Sec. 
V. 

Integrating out ./V pieces of fcrmions generates an ad- 
ditional dynamic quadratic term for the gauge field: 



d 3 k 



(2tt)< 



a^(-k)Il^(k)a v (k) 



To one-loop order: 



(5) 



(6) 



Adding the non-local gauge fixing term and combining 
the original Maxwell term in Eqn.l with the generated 
dynamic gauge field terms in Eq, 5, the propagator of the 
gauge field is [1] : 



D 



7,2 I Ne 2 i 

ft, ~r -^g rh 



(7) 



In general, the Ultra- Violet ( UV ) k — > oo behavior 
of this propagator is controlled by the original Maxwell 
term ~ 1/fc 2 , the Infra- Red (IR) one k — > is dictated 
by the generated term ~ 1/fc. In this paper, we consider 
large N limit, such that Ne 2 ^> k, so we neglect the 
original Maxwell term. The above propagator becomes: 



16 1 



Ne 2 k 



(8) 



In Landau gauge a = 0, the one-loop fermion self- 
energy Feymann diagram of the Green function Eqn.2 is 
given by: 
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Fig 1: The fermion self-energy diagram in Landau gauge 

The corresponding expression is 



E(fc) 



.16 



_fq_ 7A(fc-g) A l f _ £^ 



N J (2tt) 3 ,A1 {k - qf 



(9) 



By extracting UV divergence, we find the anomalous 
exponent rj — ^ If we choose different gauge fix- 
ing parameter a, we would get r\ = jj^z — j^pr- V is 
clearly gauge-dependent quantity and its physical mean- 
ing in any Lorentz covariant gauges which includes Lan- 
dau gauge is not evident. In the following section, we will 
calculate -q in two Lorentz non-covariant gauges: tempo- 
ral gauge and Coulomb gauge. 



III. THE CALCULATION IN TEMPORAL 
GAUGE 

As stated in the introduction, the main focus of this 
paper is the gauge-invariant Green function Eqn.4. In 
temporal gauge, the equal-space gauge invariant Green 
function is the same as the gauge dependent one [12]. 
This can be seen easily as follows. 

At equal space, Eqn.4 becomes: 



G inv {0,r) =< ^(x,0)e 



'■ /J" ao(C) d fo 



$(2,t)> (10) 



In temporal gauge a — [21], the Dirac string drops 
out, the equal-space gauge-invariant Green function co- 
incides with the equal-space conventional one: 



G mv (0,r) =< V(z,0)V(z,t) >= G(0,t) 



(11) 



The strategy is to calculate the conventional Green 
function in temporal gauge clq = and then see what 
we can say about the gauge invariant Green function. 
As shown in the last section, in Lorentz covariant gauge 
Eqn.8, there is no IR divergences to worry about, wc 
need only to regularize the UV divergence. However, in 
temporal gauge and Coulomb gauge which break Lorentz 
invariance, we run into both UV and IR divergences as 
expected. Care is needed to regularize these plaguey IR 
divergences in a physical way. 



With the notation = [k, fco) ( namely ko is along z 
direction ), in ao — gauge, we can invert Eqn.5 to find 
the propagator: 



Dij(K) 



16 1 

Ne^K 



kj k 

4/) 



(12) 



where K 2 = k% 



{kf. 



The one-loop fermion self-energy Feymann diagram is 




Fig 2: The fermion self-energy diagram in temporal gauge 

The corresponding expression is 



Y,{K) = -i 



16 



d 3 Q 1»{K-Q), 



1 



- rx i liQj , 



N J (2tt)3 n (K - Q) 



(13) 



By using standard 7 matrices algebra and suppressing 
the prefactor — ij^, we can simplify the above equation 
to: 



- 7o J 



+ 7-, 



f dH 
J (2tt 



d 3 Q (fc-g)o , cf. 
d 3 Q q 2 (k + q)i - 2q- kqi 



f (K - QYQql 



(14) 



In the following, we will try to evaluate Yi{K ) by two 
different cut-offs. 



A. Cut-off in momentum and energy Q < A 

We choose the external momentum K to be along z 
axis, namely = (0,0,1), then in spherical coordi- 
nate system, the components of integral vector Q are: 
qo = Q cos 9, qi = Q sin 9 cos (f>, q^ — Q sin 9 sin <j>, the sec- 
ond term in Eqn.14 vanishes due to the <j> integral. The 
logarithmic divergence of the first term becomes: 



7o fc o 
(2tt)2 



J dx{ 



1 - 2.t 2 )(1 + .x- 2 )logA 



(15) 



where x = — cos 9. 

The integral can be rewritten as: 



1 dx 
.1 x 2 
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y 



(16) 



As expected, we run into IR divergence at x — [22] 
which is in the middle point of the contour integral on 
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the real axis from —1 to 1. Fortunately, by physical pre- 
scription, we can avoid the IR singularity at x — by 
deforming the contour as shown in Fig. 3 



IV. CALCULATION IN COULOMB GAUGE 



-l o 
Fig 3: The contour path C to bypass the IR singularity in tem- 
poral gauge 



The divergent part in Eqn.16 becomes: 
1 dx f dz f 1 dx 

2 



r 1 dx _ r dz _ r 1 

Putting back the prcfactor 
swer: 



{x + ief 



(17) 



ijf, we get the final an- 

-^(- y )lo g A^ 7ofco ^logA (18) 

We can identify the anomalous dimension as [23] 
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We expect this is the correct anomalous dimension of 
the gauge-invariant Green function Eqn.4. Note that it 
has the same sign as that calculated in Landau gauge. 



In this section, we will first evaluate the conventional 
Green function Eqn.2 in Coulomb gauge diai — and 
speculate how to calculate the gauge invariant one Eqn.4 
in this gauge. In FQH transitions as discussed in Ref. [8], 
because the Coulomb interaction between the fermions 
breaks Lorentz invariance, Lorentz invariant gauges are 
not convenient, Coulomb gauge has to be employed to 
perform RG calculations. As in temporal gauge, we also 
run into IR divergence. However, in contrast to the 
Green function in temporal gauge, the Green function in 
Coulomb gauge does not correspond to any gauge invari- 
ant quantity, the IR divergence may not be regulariable. 
This is indeed the case as stressed in Ref. [8] in FQH tran- 
sitions contexts and as demonstrated in the following in 
QED3 context. 

With the notations = (k,k ), in Coulomb gauge, 
the propagators are 



< a (-K)a (K) > = J 

< ai (-K) aj (K) > = (S t3 ^)1 



(22) 



B. Cut-off in momentum space q < A only 



The two fermion self-energy diagrams corresponding to 
the two propagators are 



In last subsection, we choose a cut-off A in Q, in this 
subsection, we introduce a alternate cut-off A only in 
momentum space q, but integrate the frequency qo freely 
from — oo to oo. The motivation to choose different cut- 
off ( or different renormalization scheme ) is to check if 
the exponent Eqn.19 is universal. 

2x ? and extracting the logarithmic 



By setting f - , 
UV divergences, we find 



£(fc) = 



7o fc o 
4tt 2 

4^" 



dx 



1 



(1+x 4 ) 5 



36x 4 



(l+x 2 ) 4 

(1-.T 2 ) 5 1 



r 2\4 



2.r 



2x 2 
log A 



• log A 



(20) 



It can be shown that the two coefficients in front of 
7 fco and 7^ are exactly the same as expected from 
Lorentz invariance. Regularizing the integral as depicted 
in Fig. 3, we find that both are equal to: 



dx , 
1 (1 



(l-x 2 ) 5 1 



x 2 ) 4 2x 2 



16 

y 



(21) 



Putting back the prefactor — ijj- , we get exactly the 
same anomalous dimension in Eqn. 19. This demon- 
strates that it is indeed a universal constant independent 
of renormalization scheme. 



To Vo 

(a) 



Fig 4: The fermion self-energy diagrams in Coulomb gauge 

The expression for Fig. 4a is: 



.16 f d 3 Q ^(K-Q)^ Q 
X a (K) = -i— I — 370 (K _ n , 2 7ot? (23) 



N J (2tt) 3 ,u {K-QY 'V 



The expression for Fig. 4b is: 

Sfc(A-) = -i 



16 f d 3 Q 7m(^-Q)m 
Nj (2tt)3 71 (K-Q) 2 q 2 ' Q 



QiQj x 1 

-TjOfcj - — ): 



(24) 



In the following, we will try to evaluate T,(K) = 
T, a (K) + Sf,(fc) by two different cut-offs. 
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A. Cut-off in momentum and energy Q < A 

We choose the external momentum K along z axis. 
Adding the logarithmic divergences from Eqn. 23 and 
24, we find: 



7o fc o 
(2n) 2 



1 (1 - 2sV , . 
dx—^ ^ — log A 

.i 1 - x z 



(25) 



Again, we run into IR divergences at x = ±1 which 
are at the two end points of the contour integral on the 
real axis from — 1 to 1 . Unfortunately, from physical per- 
scription, we are unable to avoid the IR singularities at 
the two end points x = ±1 by deforming the contour. 
Therefore, we are unable to identify the anomalous ex- 
ponent. 



B. Cut-off in momentum space q < A only 

In this subsection, we introduce cut-off A only in mo- 
mentum space q, but integrate the frequency q n freely 
from — og to oo. 

Extracting the logarithmic UV divergences, we find: 



7ofco 



l-6x 2 + x 4 



2 (/ dx 7T— 3T7T— ^-«) l0 S A 



+ ^(/ dx 



_! {\-X 2 ){l+X 2 f 

2x 2 



IT 2 V_! (1 - X 2 ){1 + X 2 ) 2 3 



)logA (26) 



It can be shown that the two coefficients in front of 
7ofco and jiki are exactly the same as expected from 
Lorentz invariance. However, in contrast to temporal 
gauge, the two coefficients differ from that in Eqn. 25. 
This causes no disturbance, because the Green function 
Eqn. 2 in Coulomb gauge does not correspond to any 
physical quantities. As shown in the context of FQH 
transitions in Rcf. [8], this IR divergence disappears in 
physical quantities like /? function and critical exponents. 
By using the same method developed in Ref. [8], it is 
straightforward to show that the same phenomenum hap- 
pens here in the context of QED3 [25]. 

In order to calculate the gauge-invariant Green func- 
tion Eqn. 4 in Coulomb gauge, we could extend the 
method developed in the following section for Lorentz 
covariant gauge to Coulomb gauge. We expect the IR di- 
vergencies in Eqns. 25, 26 to be cancelled by the inserted 
Dirac string in the gauge-invariant Green function. 



V. LORENTZ COVARIANT CALCULATION 

In this section, we will calculate the gauge invariant 
Green function directly in Lorentz covariant gauge Eqn. 8 
without resorting to the gauge dependent Green function. 



We will also compare with the result achieved in tempo- 
ral gauge. 

The inserted Dirac string in Eqn. 4 can be written as: 
r 'a M (0^= / ' a»(x)j;(x)d d x (27) 
where the source current is: 

f^)=J c dr^S(x,-^(r)) (28) 

where r parameterizes the integral path C from x\ to x 2 . 

By combining the source current with the fermion 
current j^{x) = r ip(x)j^ip(x) to form the total current 
j^(x) = jp(x) + jp(x), we can write the gauge invariant 
Schwinger Green function Eqn. 4 as: 

G inv (x u x 2 ) = ^ J VipV^( Xl )^(x 2 )e- J ' d " x ^a^ 
■ J Va^ I ddx ( s { a vi-' le ^( x ) a »( x )) (29) 

Where Z is the partition function of QED3: 

Z = J V^Va^e- f ddx( ^ {d »- iea ^ +SM) (30) 

Because G mv (x,y) is a gauge-invariant quantity, we 
can integrate out the gauge field a M in any fixed gauge 
without affecting the final result. For simplicity, we 
choose Landau gauge a = in Eqn. 8. Extension to any 
gauge fixing parameter is straightforward and results will 
be given at the end of the section. Integrating out a M in 
Landau gauge leads to 



G mv (x u x 2 ) = ^ I VrVrri.r.i)rU,U 



J dxdx'jl(x)D liV (x-x')ji(x') 



After we expand the total current in terms of the 
fermion current and the source current, the last integrand 
becomes 



j ll {x)D lxv {x - x')j v {x') 
+ 3'l{x)D^{x-x')jl{x') 
+ 2j ll {x)D llv (x-x')f v {x') 



(32) 



The first term in Eqn. 32 is just the conventional long- 
range four-fermion interaction mediated by the gauge 
field, it leads to the anomalous exponent in Landau gauge 
given in section II: 



m = 



16 1 

JV3^ 



(33) 



The second term is equal to: 



5 



16 f d 3 kk 2 {x 2 ~x x ) 2 - {k ■ (x 2 - Xl )) 2 



N(2ir) 3 J k 3 (k- (x 2 - xi)) 2 

(1 — cos/c • (x 2 — Xi j) 



16 f A dk f 1 , 1 - x 2 
dx- 



(2tt) 2 N J k J_j_ x 



(1 — cos krx) 



TV 



= — -(-27-21ogAr + cosAr+-Ar) (34) 

7T N 2 

Where r = \x 2 — xi\ and 7 is the Euler constant. 

The linear divergence is a arti-fact of the momentum 
cut-off and should be ignored in Lorentz invariant reg- 
ularization. The cosAr is a fast oscillating factor and 
averages to zero [24] . The logarithmic factor leads to the 
anomalous exponent: 



V2 



16 1 



(35) 



which has the same sign as , but three times larger than 
We can write the third term as 

-e 2 / dxtp(x)j^ip(x) / dx' v D fl „(x - x') (36) 



Eqn. 36 is essentially quadratic in the fermions. Com- 
bining it with the free fermion action leads to: 

S= f dxip(x)(-y^d^ + e 2 7 M f dx' u D^ u (x - x'))ip(x) 

J J X\ 

(37) 

In principal, the propagator of fermion < 
Tp{xi)Tp{x 2 ) > can be calculated by inverting the 
quadratic form in the above equation, but it is not easy 
to carry out in practice. Instead we can construct pertur- 
bative expansion in real space by the following Fcymann 
diagrams in Fig. 5 



X 2 Xj X 2 

(a) 



(b) 



x ' y' 



x y 
(c) 



Being just quadratic, there is no loops in the above 
Feymann series. But there may still be potential diver- 
gences. 

The explicit expression for Fig. 5b is 



Fig 5: The perturbative expansion series of Eq.36 

The corresponding expression is: 
G{xi,x 2 ) = G (xi,x 2 ) 
- e 2 f dxG {x!,x) ( dx^-f^D^xix - x')G (x,x 2 ) 

J J X\ 

/f-X-2 
dxG (xi,x) / dx' x ^^D^ x {x - x') 

■ I dyG {x,y) f dy'^ v D va {y-y')G {y,x 2 ) + --- (38) 



d 3 qi d 3 q 2 e~ iqiX - e~ iq2X 



(2tt)3 (2tt) 3 -i{ qi -q 2 )-x 
Go(qih^x^D^ ly (q 1 -q 2 )G (q 2 ) 



(39) 



where x 2 — x\ = x. 

Suppressing the prefactor ^, we can write Eqn. 39 as 
the sum of two parts F(x) = Fi(x) + F 2 {x) with 



/d 3 k e ikx f d 3 t 



d 3 q , 7^ fc M 



) 3 ' { q(q-k) 2 



and 



+ 



+ 



-2i 



-2k ■ x-j^k^ + fc 2 7 P a: p 

q ■ xq(q — k) 2 
k ■ - k ■ q-y^Xf, 



q ■ xq{q — k) 2 



(40) 



d 3 k 



(2ir) 3 k 2 



J (27T 



d 3 q k 2 ^^ + g 2 7 AJ fc M 



(k + qfq 2 



(41) 



The logarithmic divergence in the first term in Eqn. 40 
cancels exactly that of Eqn. 41. The second term in 
Eqn. 40 is convergent. So the only possible divergence 
may come from the third term in Eqn. 40: 



/ 



d 3 q k ■ ary^fc - k ■ 97^ 
(27r) 3 q ■ xq(q — k) 2 



(42) 



It is easy to find the UV divergences in the first term 
and second term in Eqn. 42 exactly cancel each other as 
follows: 



x ' ky^x^ . x ky^x^i 
2ir 2 x 2 g + 2ir 2 x 2 



logA = 



(43) 



We conclude that there is no extra divergence from 
F(x) in Eqn. 39. To order 1/N, we have 



% = 

In all, the final anomalous exponent is 
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v = vi+m + V3 



3tt 2 N 



(44) 



(45) 



This is exactly the same as that calculated in the tem- 
poral gauge. 

It is straightforward to repeat the above calculation 
by using different gauge fixing parameter a in Eqn. 8, 
then rji, r] 2 , 773 all depend on a separately in the following 

,,.„,.. „ 16 1 8a „ 16 1 8a 16a 

way. fix - w 3^ - w^,V2 - - - jf^- 

It is easy to see that a disappears in the final answer 

V = Vl+V2+V3 = 
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There are many (1/iV) 2 corrections. 771 has the stan- 
dard correction, 772 has the correction from the propa- 
gator in Eqn.8, 773 has the correction from Fig. 5c and 
the new anomalous exponent 774 has the correction from 
Fig. 6 which is the lowest order mixing diagram between 
the first term and the third term in Eqn. 32. 




Fig 6: The lowest order mixing between the first and the third 
term 

Although evaluating all these (1/A) 2 corrections is be- 
yond the scope of this paper, we have established firmly 
our results Eqn. 45 to order 1/A. 

VI. CONCLUSIONS 

In this paper, we calculated the gauge-invariant 
fcrmion Green function in 2 + 1 dimensional QED by 
different methods. The calculations in temporal gauge 
with different cut-offs in section III and Lorentz covari- 
ant calculation with different gauge fixing parameters a 
in section V all lead to the same answer. These facts 
strongly suggest Eqn. 19 is the correct exponent in the 
leading order of 1/N. These methods have been applied 
to many very different physical systems [20] . 

Finally, we would like to comment on the possible ap- 
plication of the results of this paper on high temperature 
superconductors. As explicitly demonstrated in Refs. 
[14,15], in the vortex plasma regime in the neighborhood 
of the finite temperature phase transition boundary, the 
relevant gauge- invariant fermion Green function is that of 
fcrmions moving in classical random magnetic field gener- 
ated by the vortex plasma. Obviously, this static gauge- 
invariant Green function is expected to show rather dif- 
ferent behaviors than the dynamic one calculated in this 
paper. For technical and simplicity reason, the dynamic 
gauge-invariant one was used in Ref. [14] to present the 
Energy Distribution Curve (EDC) and Momentum Dis- 
tribution Curve (MDC) of ARPES data. We used the 
anomalous exponent calculated by RW which differs from 
our result Eqn. 19 both in sign and in magnitude. How- 
ever, if correct anomalous exponent is used in the calcu- 
lation in Ref. [14] is irrelevant anyway, because as cau- 
tioned at the end of Ref. [14], the static gauge- invariant 
Green function should be used to compare with experi- 
mental ARPES data [20]. 
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Note added in proof. After we achieved the re- 
sults of this paper, we got to know the preprint cond- 
mat/0112202 by Khvcshchenko. He used very different 
methods and got the same anomalous exponent as mine 
Eqn. 19 upto a factor of 2. I suspect that the factor of 2 
discrepancy may be due to he used 4 component spinors 
in contrast to the two component spinors used in this pa- 
per. Four component spinor representation has the ad- 
vantage to investigate dynamical parity-conserving mass 
generation as discussed in Refs. [1,3]. 
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